Soliton ratchets 
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The mechanism underlying the soliton ratchet, both in absence and in presence of noise, is investi- 
gated. We show the existence of an asymmetric internal mode on the soliton profile which couples, 
trough the damping in the system, to the soliton translational mode. Efi^ective soliton transport is 
achieved when the internal mode and the external force are phase locked. We use as working model 
a generalized double sine-Gordon equation. The phenomenon is expected to be valid for generic 
soliton systems. PACS numbers: 05.45.Y, 05.45.-a, 05.60. C, 63.20.Pw 



One of the phenomena which is presently attracting 
interest both in physics and in biology is the so called 
ratchet effect [|j. In simple terms, a ratchet system can 
be described as a periodically forced brownian particle 
moving in an asymmetric potential in presence of damp- 
ing and periodic driving M. The periodic forcing keeps 
the system out of equilibrium so that the thermal en- 
ergy can assist the conversion of the ac driver into ef- 
fective work (direct motion of the particle) without any 
conflict with the second law of thermodynamics. This 
phenomenon has been found in several physical 0] and 
biological El systems and is presently considered as a 
possible mechanism by which biological motors perform 
their functions g. For ode systems with damping, ad- 
ditive forcing and noise, the ratchet effect can be viewed 
as a phase locking phenomenon between the motion of 
the particle in the periodic potential and the external 
driver [|6|. Ratchet dynamics have been observed also in 
more complicated systems such as overdamped (j)'^ mod- 
els [0 , chains of coupled particles with degenerate on-site 
potentials S, long Josephson junctions with modulated 
widths 1^ , inhomogcncous parallel Josephson arrays [Q , 
3D helical models |ll| , etc. These are infinite dimensional 
systems described by continuous or discrete equations 
of soliton type, with asymmetric potentials, damping, 
and periodic forcing, in which the ratchet phenomenon 
manifests as unidirectional motion of the soliton (soli- 
ton ratchet). For overdamped systems one can reduce 
the soliton ratchet to the usual single particle ratchet by 
using a collective coordinate approach for the center of 
mass of the soliton . For underdamped or moderately 
damped systems, however, this could be unappropriate, 
since the radiation field present in the system can play 
an important role for the generation of the phenomenon. 

The aim of this Letter is to investigate the mecha- 
nism underlying soliton ratchets both in absence and in 
the presence of noise. To this end we use an asymmet- 
ric double sine-Gordon equation as a working model for 
studying the effect (the phenomena, however, will not 
depend on the particular model used). We show that 
the asymmetry of the potential induce a spatially asym- 
metric internal mode on the soliton profile which can be 
excited by the periodic force. In presence of damping. 



this mode can exchange energy with the translational 
mode so that the soliton can have a net motion under 
the action of the ac force. In this mechanism, the damp- 
ing plays the role of coupling between the internal mode 
and the translational mode of the soliton. We find that, 
for fixed amplitude and frequency of the ac force, there 
is an optimal value of the damping for which the trans- 
port (i.e. the velocity achieved by the soliton) becomes 
maximal. In this case the frequency of the internal mode 
and the one of the external driver, perfectly match (phase 
locking). A similar resonant behavior is also observed by 
varying the frequency of the ac force, keeping fixed the 
other parameters of the system. At very high damping 
and fixed amplitude of the forcing, the internal mode os- 
cillation becomes very small, and the transport due to 
soliton-ratchct is strongly reduced. At low damping and 
higher forcing we find, quite surprisingly, that current re- 
versals can occur. Finally, we show that soliton ratchets 
can survive the presence of noise in the system. 

We start by introducing the following generalized dou- 
ble sine-Gordon equation 



= - sin(0) - A sin(20 + 0) = - 



dU{(t)) 



(1) 



with the potential C/((/)) = C - cos((/)) - | cos(20 + 9). 
Here A is the asymmetry parameter, ^ is a fixed phase 
and C a constant which fix the zero of the potential. A 
discrete version of this equation was introduced in Ref. 
Q in terms of a chain of elastically coupled double pen- 
dula assembled by a gear of ratio 1/2 with a phase an- 
gle 9 between them. For A = Eq. (|l|) gives the well 
known sine-Gordon equation (SGE) with exact soliton 
solutions, while for A 7^ and 9 — (mod it) it re- 
duces to the proper double sine-Gordon equation (note 
that in both cases the potential is periodic and symmet- 
ric in 4>). In the following we are interested in the case 
9^0 (mod tt) for which the periodic potential becomes 
asymmetric. We shall refer to this case as the asymmet- 
ric double sine-Gordon equation (ADSGE). In particu- 
lar, we fix 6* = 7r/2 in Eq. (1^) in order to have maxi- 
mal asymmetry, and choose C = cos (^0) — ^/2 sin (2(/)o), 
with 00 = arcsin[(l - A)/4A] -I- 2n7r and A = Vl + SA^, 
to have the zero of the potential in correspondence with 



its absolute minima <f>o. We also assume, for simplic- 
ity, As [— 1 , 1] to avoid relative minima appearing in 
the potential. Besides the mentioned mechanical model, 
Eq. (|l|) is also linked to another interesting physical sys- 
tem i.e. a one dimensional array of inductively coupled 
SQUIDs, each consisting of a loop of a Josephson junction 
in parallel with a serie of two identical Josephson junc- 
tions. The single element of this array was studied in Ref. 
p3[ in which it was shown that, due to the ratchet effect, 
the system can rectify periodic signals. In analogy with 
the single particle ratchet, it is reasonable to introduce 
in the distributed model, periodic forcing, damping, and 
noise, this leading to the following perturbed ADSGE 

(f>tt - <t>xx + sin(0) -I- A cos(2(^) = (2) 

— a4>t + esin(a'i -I- Oq) + n{x, t). 

Here a denotes the damping constant, n{x,t) is white 
noise with autocorrelation 



< n{x, t)n{x',t') >= D5{x - x')5{t - t'), 



(3) 



and e, w, ^o are, respectively, the amplitude the frequency 
and the phase of the driver. Travelling wave solutions of 
the unperturbed ADSGE, i.e. solutions which depend 
on ^ = (a: - Vt)/yjl - V^ , (note that Eq. (0) is Lorentz 
invariant), can be found by substituting </> = (/)(^) in the 
left hand side of Eq. (pi) and equating it to zero. After 
one integration in ^ we obtain 



d<t) 



00 ^2{E + U) 



e-Co, 



(4) 



which gives, after the inversion of the integral at _E = 
(top of the reversed potential), the 27r-kink (antikink) 
solutions as 



^1 



2tan"^{- 



sign{\)AB 



A — \ — ri sinh 



±l^AB/\\\ 



} (5) 



where 77 = 2Av^2(H-yl), and B = ^2(4A2 - I + A) (the 
plus and minus signs refer to the kink and antikink so- 
lutions, respectively). Note that in the limit A — > Eq. 
(H) reduces to the well known soliton solution of SGE. To 
investigate the existence of internal modes in the system, 
we linearize the ADSGE around the solution in Eq. (|^), 
i.e. we look for solutions of the form = (jyf. + ip with 



V'(a;,i) = exp(jwt)/(x), f{x) « 1. 

This leads to the following eigenvalue problem on the 
whole line 



fxx + (^' - cos[0|] + 2Asin(20|)) / = 0, 



(6) 



with fx{±oo) — 0, which can be easily solved by numer- 
ical methods for any finite length of the system. In Fig. 
|l| we report the numerical spectrum of Eq. (H) as a func- 
tion of A. We see that except for the sine-Gordon limit 



(A = 0), there is an internal mode frequency 17^ below the 
spectrum of the phonon band (the zero mode, existing for 
all values of A, is not plotted for graphical convenience). 
In the inset of the figure the shape of the internal mode 
is also reported, from which we see that the asymmetry 
of the potential induces a spatial asymmetry in f{x). In 
the presence of a periodic force this internal mode can be 
easily excited. 
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FIG. 1. Small oscillation spectrum versus A. Above the 
solid line are the phonon's modes. The starts joined by the 
dashed line refers to the internal mode. The inset shows the 
spatial profile of the internal mode at A = 0.5. 

To understand the role of the various elements of the 
problem, (i.e. asymmetry of the potential, internal mode, 
damping, forcing, and noise) it is better to consider first 
the zero noise case (deterministic soliton ratchet). By 
viewing the soliton as a string lying on the potential sur- 
face S'( [/, 0, a;) and connecting adjacent minima, the fol- 
lowing picture of the phenomenon can be given. If the 
potential is asymmetric (in 0) the transition from the top 
of the potential to one minimum and from the top of the 
potential to the other minimum, is also asymmetric (it 
will be more rapid for the part of the string lying on the 
region where the potential is more stiff). Thus, the po- 
tential asymmetry in </> induces an asymmetry in space 
which can be seen both in the 27r-kink profile and in the 
internal mode. In presence of an ac force, but in absence 
of damping, this asymmetry will not produce transport, 
i.e. the string will slide, without any friction, back and 
forth on the potential profile along the x— direction. The 
presence of damping, however, introduces friction in this 
sliding, and the part of the string moving on the stiff 
part of the potential profile dissipate more that the other. 
This asymmetry in the dissipation produces net motion 
for the soliton (the string moves in the direction in which 
it approaches the potential minimum more smoothly). 
We can say that the effect of the damping is to couple 
the internal mode to the translational mode. The mech- 
anism underlying deterministic soliton ratchet can then 
be described as follows: the ac force pumps energy in the 
internal mode which is converted into net dc motion by 



the coupling with the zero mode induced by the damp- 
ing. ^From this picture one can easily predict that in 
absence of the internal mode, or in absence of damping, 
no soliton ratchet can exist. Moreover, one expects that 
the maximal effect in transport, is observed when the in- 
ternal oscillation and the external force are synchronized 
(phase-locked) . 

In order to confirm this picture, we have performed 
direct numerical integrations of the ADSGE for different 
values of the system parameters. First, we have checked 
that in the SGE limit, i.e. when A = 0, the ratchet dy- 
namics does not exists. This agrees with the fact that 
in this case there is no asymmetry in the potential and 
no internal mode. At this point we remark that the dy- 
namics of a SG kink subject to a periodic force was also 
investigated in Ref. |14|, in which it was shown that in 
absence of damping the kink can acquire a finite veloc- 
ity depending on the initial phase of the ac force. Net 
motion of a SGE soliton was shown to be possible also 
in presence of a small damping, if the ac force excites 
a phonon mode which exchange energy with the soliton 
[ p^ . These cases, however, should not be confused with 
soliton ratchets since they strongly depend on initial con- 
ditions (if one average on initial conditions the transport 
disappears). Moreover, in contrast with soliton ratchets, 
these effects exist only at zero or at very low damping. 
Second, we have checked that for A ^ (asymmetric po- 
tential) but in absence of damping, soliton ratchets also 
do not exist (for brevity we will not expand on these cases 
here). From this analysis we conclude that, in analogy 
with the deterministic single particle ratchets, the asym- 
metry of the potential, the damping and, obviously the 
periodic forcing, are crucial ingredient for soliton ratch- 
ets. In Fig. y we show a prospectic view of the soliton 
ratchet dynamics as obtained from numerical integration 
of Eq. (||). We remark that the direction of the motion 
is fixed by the asymmetry of the potential and can be 
inverted by changing the sign of A. We see that, except 
for the soliton profile which is wobbling, no phonons are 
present in the system. This confirms the relevance of 
the internal mode in the phenomenon. We also find that 
for some parameter values, the net motion is more effec- 
tive. To investigate the dependence of the phenomenon 
on parameters, we have performed numerical simulations 
of Eq. (0) both by fixing all parameters and changing w, 
and by fixing all parameters and changing a. In Fig.0 we 
show the average velocity (computed by using an inte- 
gration time t — 1000) of the soliton center of mass ver- 
sus the frequency of the external driver, for two different 
values of the amplitude of the ac force (the low damping 
part of the curves was not computed due to the longer 
integration times required in this case). From this figure 
we see that {V) has a maximum at uj ^ 1, (i.e uj = 1.04 
and w = 1 for e = 0.4, e = 0.6, respectively), this being 
very close to the internal mode frequency fli ~ 1.0562 
(the discrepancy is within the numerical accuracy of our 
numerical scheme). By increasing the amplitude of the 
forcing, the dynamics gets more complicated (breather- 



like excitations can appear) and the resonance peak in 
frequency more pronounced. 




FIG. 2. Time evolution of the x-derivative of the kink pro- 
file while executing ratchet dynamics. The parameter values 
are A = 0.5, Oo = 7r/2, a = 0.6, e = 0.8, ui = 0.8, and D = Q. 
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FIG. 3. Mean velocity of the kink center of mass versus the 
frequency of the ac force for e = 0.4 (solid line) and e = 0.6 
(dashed line). The other parameter are fixed as in Fig.(Q) 
except for the asymmetry parameter which is A = —0.5. 

A similar resonant behavior is expected to exists also 
as function of the damping. When the damping is very 
high, indeed, the internal mode is almost suppressed by 
the damping, while when it is very low the coupling be- 
tween the internal mode and the translational mode is 
very small, both cases giving minimal transport. In be- 
tween these extremes, a value of the damping which allow 
the internal mode to synchronize with the external driver 
and optimize transport, should then exist. This is what 
we observe in Fig. where the velocity vs damping is 
reported as a continuous curve. Although the ratchet ve- 
locity can be slightly increased by optimizing parameters, 
it remains small in comparison to the velocities achieved 
under the action of dc forces. This can be due to the 



fact that the internal and the translational mode seems 
to couple only to second order in a perturbation scheme. 
To increase the velocity one could increase the ampli- 
tude of the ac force, but this is limited by the threshold 
above which the system becomes chaotic. In practical 
applications, however, a small drift velocity is sufficient 
to manipulate solitons under the action of external ac 
forces. 
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FIG. 4. Mean velocity of the center of mass of the kink 
versus q, for parameter values A = —0.5, ^o = 7r/2, uj — 0.4, 
e = 0.5. The continuous curve refers to the case D = 
while the stars refer to the nondeterministic case D = 0.01. 
The continuous and dotted curves in the inset refer to two 
soliton (px profiles at D = 0, taken in the co-moving frame 
and separated in time by one period of the driver. 

To show that the internal mode is phase locked with 
of external driver we have plotted in the inset of this fig- 
ure the (f)x profiles in the co-moving frame (i.e. the drift 
motion was subtracted) at two fixed times ti = 211.6 
(solid line) and ^2 = 227.3 (circles) separated by one 
period T — 2Tr/uj = 15.7 of the driver. We see that 
the profiles overlap each other, i.e. the oscillation on 
the kink profile is perfectly synchronized with the exter- 
nal driver (phase locking). /From Fig. ^ we see, quite 
surprisingly, that at low damping current reversals can 
occur (note that the average velocity becomes negative 
for a less than acr ^ 0.24). We find that the value of 
acr increases as the amplitude of the driver is increased. 
The occurrence of this phenomenon, which resemble the 
one observed in single particle ratchets at low damping 
P,p^ , seems to be related more to the phonon-soliton in- 
teraction than to the internal mode mechanism described 
above (at low dampings a complicate transferral of en- 
ergy between phonons, internal mode, and translational 
mode, can arise) . To understand this phenomenon, how- 
ever, a more detailed study is required. 

Finally, we have investigated the effect of the noise 
on deterministic soliton ratchets. A preliminary analysis 
shows that for low noise amplitudes the effect of the noise 
on the phenomenon is minimal, in the sense that the dy- 



namics gets dressed by the noise, but after averaging on 
the noise, almost the same soliton mean velocity {V) is 
obtained. This is shown by the stars in Fig. ^ which rep- 
resent the numerical values of (V) calculated in presence 
of a noise of amplitude D — 0.01. The fact that soliton 
ratchets can survive the presence of weak amplitude noise 
can be understood as consequence of the structural sta- 
bility of phase locking phenomena against "small" fluc- 
tuations. This indicates that the phenomenon can exist 
also in real systems such as one dimensional arrays of 
inductively coupled SQUIDs. The mechanism of soliton 
ratchets discussed in this Letter is expected to be valid 
also for other soliton systems. 
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